Quantum-Fluctuation-Driven Coherent Spin Dynamics in Small Condensates 
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We have studied quantum spin dynamics of small condensates of cold sodium atoms. For a 
condensate initially prepared in a mean field ground state, we show that coherent spin dynamics are 
purely driven by quantum fluctuations of collective spin coordinates and can be tuned by quadratic 
Zeeman coupling and magnetization. These dynamics in small condensates can be probed in a high- 
finesse optical cavity where temporal behaviors of excitation spectra of a coupled condensate-photon 
system reveal the time evolution of populations of atoms at different hyperfine spin states. 



Recently, single-atom detection in optical cavities has 
been realized in experiments by having atoms and cav- 
ity photons in a strongly coupling regime [H, This re- 
markable achievement has been applied to study optically 
transported atoms in cavities furthermore the cou- 
pling between a small Bose-Einstein condensate (BEC) 
and cavity photons and resultant collective excitations 
have also been successfully investigated The sensitiv- 
ity that a cavity-based atom detector has, together with 
a translating optical lattice which can effectively trans- 
port ultra cold atoms from a magnetic-optical trap to 
a cavity make it possible to study the physics of small 
BECs. Especially, this potentially opens the door to 
explore coherent dynamics of ultra-cold atoms in rela- 
tively small condensates. The physics of BECs of small 
numbers of atoms can be qualitatively different from 
the physics of big condensates and represents a new do- 
main of cold-atom research. In small condensates, vari- 
ous intrinsic beyond-mean-field dynamics can be relevant 
within an experimentally accessible time scale. These 
new physical phenomena however have been quite diffi- 
cult to study using the standard absorption-imaging ap- 
proach to cold atoms because of relatively fewer atoms 
are involved in small condensates. Cavity electrodynam- 
ics in a strong coupling regime and high sensitivities to 
intra-cavity atoms on the other hand are ideal for in- 
vestigating small condensates where the beyond-mean- 
field dynamics are mostly visible. In this letter, we fo- 
cus on the basic concepts of beyond-mean-field coherent 
spin dynamics in BECs with typically a few tens to a 
few hundreds of atoms and detailed analysis of detect- 
ing these fascinating properties of small condensates in 
optical cavities with high-finesse. Research on this sub- 
ject could substantially advance our understanding of the 
nature of quantum- fluctuation dynamics^, in this par- 
ticular case, dynamics purely driven by fluctuations with 
wavelengths of the size of condensates. Secondly, results 
obtained can help to better recognize limitations of pre- 
cise measurements of various interaction constants based 
on mean- field coherent dynamics Q. Thirdly, our results 
should shed some light on the feasibility of investigating 
fluctuation dynamics of small condensates using optical 



cavities and also pave the way for future studies of dy- 
namics of coupled small condensates. 

To understand spin dynamics of a small condensate, we 
first study the evolution of a condensate of N hyperfine 
spin-one sodium atoms which is initially prepared in a 
mean field ground state. 
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Here n is a unit director and three components of ■0^, 
^l^, a = x,y, z are creation operators for three spin-one 
states, \x) = (|1) - I - 1))/V2, \y) = (|1) + | - l))/iV2 
and \z) = |0) respectively. And in this representation, 
Sa = —i^aPjipJ^'^'j is the total spin operator. States in 
EqUwith n = 62 minimize the interaction energy of the 
following Hamiltonian for spin-one atoms in the presence 
of a quadratic Zeeman coupling along the z-direction. 
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Here C2 is a spin interaction constant and q is the 
quadratic Zeeman coupling 0, i, i Hi. Mean field 
ground states are stationary solutions to the multi- 
component Gross-Pitaevskii equations for spin-one atoms 
and dynamics of these initial states demonstrated below 
are therefore a beyond- mean field phenomenon. When 
deriving Eq[5] for a trapped condensate, we assume 
that spin dynamics are described by a single mode, i.e. 
lAalr, i) — \/ p{r)^a{t)', for a small condensate of less 
than one thousand weakly interacting atoms, this ap- 
proximation is always valid. C2 is typically a few nano 
kelvin for sodium atoms; q — {fiBBy / {AAhf) and the 
hyperfine splitting is Ahf = (2Tr)1.77GHz [fiB is the 
Bohr magneton and h is set to be unity). 

To illustrate the nature of non-mean-field dynamics 
and crucial role played by quantum fluctuations, we ex- 
pand the full Hamiltonian in EqO around a mean field 
ground state. In the lowest order expansion, we approx- 
imate ■0''' « VNgz + '4>x^x + ipl^yi and 'ipl:,y a-re much 
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less than V N; the Hamihonian then can be expressed in 
terms of the bihnear terms 
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where for a — 
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"00,) and Pa = 



ipa) are pahs of conjugate operators which 
satisfy the usual commutation relations [9a, Pp] = iSa,i3- 
Semiclassically, collective coordinates 9a, a = x,y rep- 
resent projections of -0^ or order parameter n in the xy 
plane, and Px{y) ~ Sy^x) is the spin projection along the 
?/(x)-direction. The bilinear Hamiltonian is equivalent to 
a harmonic oscillator moving along the direction of 9x.y 
with a mass ?ne// = q+ic2 ' ^ harmonic oscillator fre- 
quency w — \J q{q + 4c2) and effective spring constant 
qN; the mass at g = is induced by scattering between 
atoms. The excitation spectrum is En — {n + 1/2)lu, 
n = 0, 1, 2... . When g = 0, the Hamiltonian describes a 
particle moving in a free space. 
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FIG. 1: a) (color online) Time evolution of No{t), the atom 
population at |1,0) state for different quadratic Zeeman cou- 
pling q. Initially, all TV = 200 atoms occupy |1, 0) state which 
corresponds to a mean field ground state. The inset is tc, the 
time for the first dip in the 5 = data, as a function of the 
number of atoms A'', b) Time evolution of No{t) for different 
magnetization m (here q = 0). All initial states are again 
mean-field ground states for given m. Inset is for m = 0.98. 
In this and FigllEl C2 = {2tt)50Hz. 

In the ground state of the bilinear Hamiltonian of 
EqEl {9a) = {Pa) = and n and (S) have no projec- 
tions in the xy plane. However, quantum fluctuations of 



^aj.y-coordinates in the ground state can be estimated as 
{9a9a) = \J This is a measure of how strongly 
n fluctuates in the xy-plaxie. As expected, these quantum 
fluctuations are substantial only when q is small and are 
suppressed by a quadratic Zeeman field which effectively 
pins the order parameter along the z-direction. A direct 
calculation also shows that the amplitude of quantum 
fluctuations {9a) mf in the mean field ground state de- 
fined in Eq[T]is 1/2N. This indicates that the mean field 
ground state is a good approximation only when q ^ 4c2. 
On the other hand, as q decreases and the effective spring 
constant gets smaller, the deviation becomes more and 
more severe. When q approaches zero, quantum fluctua- 
tions 9a in the harmonic oscillator ground state become 
divergent implying that the mean field ground state is no 
longer a good approximation. 

Indeed, the the energy of mean field ground state is 
Emf = 2 "I" "^2 which is much higher than ^cu when q 
C2; such a state corresponds to a highly excited wave 
packet, because of a relatively narrow spread along Re- 
directions and consequently an enormous kinetic energy 
associated with momenta Pa- We therefore expect that 
dynamics in this limit could dramatically differ from a 
stationary solution. Since the total number of atoms N 
is equal to ipa'^^a, the population of atoms at \z) (or 
|1,0)) state A^o = (V'lV'^) 

is directly related to quantum 

fiuctuations of 0„ and Pa, 
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EqH] shows that the time evolution of NQ{t) is effectively 
driven by quantum fiuctuations in 9a and Pa ', a study of 
No(t) probes underlying quantum- fluctuation dynamics. 

For an initial state prepared in a mean field ground 
state with n — where all atoms condense in 1 1 , 0) 
state, one finds that {9"^) = and (P^) = y- The evo- 
lution of such a symmetric Gaussian wave packet subject 
to the bilinear Hamiltonian can be solved exactly using 
the standard theory for harmonic oscillators. The wave 
packet will remain to be a Gaussian one with the width 
oscillating as a function of time. Qualitatively, because 
of the symmetry, only harmonic states with even-parity 
are involved in dynamics and therefore the oscillation 
frequency is 2lu. Furthermore during oscillations, the ki- 
netic energy stored in initial wave packets is converted 
into the potential one and vice versa. Especially when 
q ^ C2, oscillations are driven by the enormous initial 
kinetic energy associated with Pa', the oscillation ampli- 
tude can be estimated by equaling the total energy Em f 
to the potential energy which leads to {9a) ~ C2/{Nq). 
A straightforward calculation yields the time dependence 
of {9l) and {P^) that leads to 
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The oscillating term in Eql5] shows the deviation from 
the stationary behavior due to quantum fluctuations in 
0Q,-coordinates. The deviation is insignificant when q is 
not too small; however when q is of the order of 02/ N, we 
expect that the non-mean held dynamics becomes very 
visible. Note that the approach outlined here neglects 
all higher order anharmonic interactions and therefore is 
only valid when the relative amplitude of fluctuations is 
small; that is when q ^ C2/N. 

When q approaches zero, the short time dynamics fol- 
lowing the bilinear Hamiltonian is equivalent to a particle 
of a mass nief / = N/Ac2 that is initially localized within 
a spread (0^) = 1/2N having a ballistic expansion with 
a typical velocity give as (u^) = 8C2/N. The time de- 
pendence of spread (0^) therefore is 1/2N + {%c\/N)t^. 
So at t '--^ -\//V/c2, the number of atoms not occupy- 
ing the initially prepared |1,0) state becomes of order 
of N . This limit was first addressed by Law et al. in the 
context of four- wave-mixing theory and also in early 
works 11 , l2| ; to describe the physics after this charac- 
teristic time scale requires analysis of full quantum dy- 
namics. This time scale however becomes quite long for 
a few million atoms which makes it difficult to observe 
quantum dynamics in large condensates. 

In the following, we are going to present our numeri- 
cal results on dynamics and focus on its dependence on 
quadratic Zeeman coupling q and magnetization m. For 
a condensate of iV = 200 atoms, we numerically integrate 
the time-dependent N-body Schrodinger equation of the 
quantum Hamiltonian in Eq|21 The time evolution of A^o 
driven by quantum fluctuations is shown in Fig[T^). As q 
increases far beyond 0.2c2, A^o oscillates as a function of 
time with frequency 2w and the amplitude of oscillations 
decreases; the damping is not visible over tens of oscilla- 
tions. When q is below 0.2c2, anharmonic effects become 
substantial and oscillations are no longer perfect; when 
q = C2/4O, oscfllations are strongly damped after a few 
cycles and revived afterwards. For q — 0, Nq drops to a 
minimum of about 0.38N when t = — 0.53\/^ /c2 and 
remains to be a constant before reviving to be 0.8A^ at 
about lOic. For sodium atoms with a typically density 
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C2 = (27r)50i7z; = 23.8tos for TV = 200 



and increases to a few seconds when iV reaches 2 x 10^. 

We have also studied the quantum dynamics of a mean 
field condensate with a finite magnetization along the z- 
direction, m = me^. States which minimize the mean 
field energy of the Hamiltonian in Eq[2] with (7 = are 
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where svcl2,j] = m, m(€ [—1, 1]) is the normalized magne- 
tization. By expanding the Hamiltonian around these 
mean field states, one obtains a harmonic oscillator 
Hamiltonian defined in terms of conjugate operators, 
^ ^(^t + and = - ^z)- The effec- 

five mass is mef f = — -. — ^. and the harmonic os- 

cillator frequency lu = 2|m|c2. States shown in Eql6]have 
a narrow width along the direction of 9^, (6*^) = 1/2A^ 
and therefore carry large conjugate momenta Pz', the cor- 
responding large kinetic energy drives a unique non-mean 
field quantum spin dynamics. The harmonic expansion 
again is only valid when m is large and fluctuations are 
weak. Simulations of the full Hamiltonian have been car- 
ried out in this case; in Fig[T]D), we show the time depen- 
dence of No{t) for different magnetization. Only when m 
is close to unity, weakly damped oscillations are observed. 
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FIG. 2: (a,b,c) (color online) Eigenfrequencies Ap as a func- 
tion of t for different detuning Ac when the relative population 
at state |1, 0), po{t) evolves. At i = 0, all atoms occupy |1, 0) 
state and N = 200. d) Ap as a function of time t for Ac — 0. 

We propose a method to probe quantum spin dynam- 
ics of a small condensate of spin-one sodium atoms using 
cavity quantum electrodynamics. For a Bose-Einstein 
condensate with N atoms coupled to a quantized field 
of a cavity, a single cavity photon can coherently inter- 
act with atoms which leads to a collective coupling of 
g-y/iV [l3|. In experiments [1, 0], atoms are transported 
into a cavity via a moving optical lattice; excitations are 
measured by individual recordings of cavity transmission 
when frequencies of an external probe light are scanned. 
Here we consider a multi-component BEC coupled to a 
single cavity mode; the eigenfrequencies of the coupled 
system uniquely depend on populations at three hyper- 
fine states. By measuring the energy spectrum of this 
coupled system, one obtains temporal behaviors of atom 
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populations at different states. 

We restrict ourselves to excitations which involve a sin- 
gle cavity photon interacting with atoms in a BEC. We 
study atomic transitions from 35*1 — > 3Pi in sodium 
atoms. The Hamiltonian consists the following terms, 

Heavily = T^Wg^glgi + hwe^elej + TiWcclcp 

i j P 

5Z dfjSjCpgt + h.c, (7) 

P i,3 

where i labels three states \F = l,mp = 0,±1) in 3Si 
orbital and j eight states \F' = 1,to_f'), \F' = 2,mF') in 
3Pi orbital. g| and el create atoms in one of 3Si and 
3Pi states respectively with corresponding frequencies 

Wg- ,Wej ■ c|, creates a photon with frequency Wc and po- 
larization p in the cavity mode. gfj{— Dfjy^hwc/2eQV) 
is the coupling strength for a transition i j driven by a 
cavity photon with polarization p, which depends on the 
dipole matrix element D^-, the effective mode volume V. 
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FIG. 3: Eigenfrequencies Ap as a function of time t driven 
by dynamics of population po{t) in the presence of quadratic 
Zeeman coupling q = 0.05c2 (or B = 95mG). Again at t = 0, 
all atoms occupy |1, 0) state and TV — 200, Ac = 0. 

For simplicity, we set the energy of 35*1 states to 
be zero, i.e. Wg^ = 0; the energy of excited states is 
Wg'^ = Wa ± A with 2A being the hyperfine splitting 
between F' — 2 and F' = 1 states. For atomic tran- 
sitions induced by left-circularly(cr"'") polarized photons, 
the selection rule is AF = 0,±1, Am^ = 1- In a cavity, 
a state with a cavity photon (labeled as lc),iV„i^ atoms 
at 3Si lljUip) states and no atoms at excited states (la- 
beled as Oj) is expressed as \ lc; Ni, Nq, N-i;Oj); it is cou- 
pled to the following states with one of atoms excited to 
3Pi states (labeled 1_f') and no cavity photons (as Oc), 
\oJNi - 1, No, 7V„i; 1f'=2), \0c; Ni,No- 1, TV-i; 1f'=i), 
|0e; Ni,Nq- 1, iV_i; 1f'=2), \0c] Ni,No, TV_i - I; l^^^i), 
|0c; iVi, A^O) -^-1 — 1; 1_F'=2)- We diagonalize the Hamil- 
tonian matrix and obtain six eigenfrequencies ujp for this 
coupled system. Three are 3Pi orbitals without mix- 
ing with 351 states, Wp — Wa ^ A; the other three de- 
pend on relative populations of atoms at each spin state. 



Pmp = Nmp/N, mp = 0,±1. The latter three eigen 
frequencies are determined by the eigen value equation, 

(Ap - Ae)(A2 - A2) - ApNglFi - ANglF^ - 0. (8) 

Here Ap = Wp — Wa, m = p+i — p_i is the normal- 
ized magnetization, and gi is the coupling between 35*^ 
|F = l,mF = 1) and 3Pi \F' = 2,tof' = 2) by cr+ 
light; Pi = (2 + to)/3 and Pj = (1 + m)/2 - po/6. Ap- 
parently eigenfrequencies Api^p2,p3 are a function of pmp 
and therefore can be used to probe the variation in po,±i 
due to coherent spin dynamics. Api^p2,p3 depend on a 
dimensionless parameter r — "■^si _ when detuning 

Ac — and as r ^ oo, Api,p2,p3 are around 0, ±-^VlVgi 
respectively. Eigenfrequency Ap2 varies from —3 A/4 to 
—A/2 when po increases from to 1; the variation am- 
plitude S = Ap(po = 1) — Ap(po = 0) reaches a satu- 
rated value A/4. For sodium atoms, A = {2^)94 AM Hz; 
cavity parameters are chosen according to Ref.Q] and 
51 = {2Tr)10MHz. In Fig. ([2]), we show the evolution 
of Ap in time for different detuning Ac when atoms are 
initially prepared at state |1,0) of 351/2. The evolution 
of Ap{t) which can be probed by a beam maps out 
population No{t) driven by underlying quantum fluctu- 
ations. In Fig. ([3]), we further show the time depen- 
dence of Ap due to oscillatory quantum spin dynamics 
for q = 0.05c2(P « 95toG), N = 200 and A^ = 0. 

In conclusion, we have illustrated the nature of co- 
herent spin dynamics driven by quantum-fluctuations in 
small condensates. The time evolution of population of 
atoms at different hyperfine spin states is shown to re- 
veal intrinsic dynamics of quantum fluctuations of or- 
der parameters and spin projections. These dynamics 
can be probed by studying eigenfrequencies of a coupled 
condensate-photon system in a high-finesse optical cav- 
ity available in laboratories. We thank Gerard Milburn, 
Junliang Song for stimulating discussions. This work is 
in part supported by NSFC, 973-Project (China), and 
NSERC (Canada), CIFAR and A. P. Sloan foundation. 
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